ASYMPTOTIC BEHAVIOR OF RANDOM HEAPS 



J. BEN HOUGH 1 

Abstract. We consider a random walk W n on the locally free group (or equivalently a signed random heap) with 
m generators subject to periodic boundary conditions. Let #T(W n ) denote the number of removable elements, 
which determines the heap's growth rate. We prove that lim n _oo - w n " < 0.32893 for m > 4. This result 
disproves a conjecture (due to Vershik, Nechaev and Bikbov [4]) that the limit tends to i as m — > oo. 



1. Introduction 

We consider a random walk W n on the locally free group hT m with to generators subject to periodic boundary 
conditions, (go, g±, ... , g m -i '■ gigj — gjgi V i =/= j ' ± 1 mod to). The random walk on this group has an elegant 
interpretation in terms of random heaps, which are defined in P] and examined more carefully in g]. Consider 
the lattice Z m x Z + , where Z m denotes the integers modulo to, and drop signed pieces with +/— charges 
uniformly over the m columns. When a piece is dropped in a column it falls as far as possible subject to the 
condition that it cannot fall past a piece in its own column or a piece in either of the neighboring columns, 
i.e. the pieces have "sticky corners" . Also, if a piece with a + charge lands directly on top of a piece with 
a — charge, the two pieces annihilate. More precisely, each piece in the random heap may be described by 
three coordinates: its height, horizontal position and sign. Let(/ij, j, Sj) denote the highest piece in column 
j, and set hj = if no such piece exists. When a new piece is to be added to the heap, its horizontal position 
k and sign s are chosen uniformly over the 2m possiblities. If columns k — 1, k and k + 1 are empty, then 
the height of the new piece is 1. If hk > max{/ifc_i, hk+i} and Sk ^ s, the new piece and (hk,k,Sk) will 
both annihilate. Otherwise, the new piece is added to the heap and its height is max{/in, hk, hk+i} + 1. 
A typical random heap is illustrated in Figure ^ 
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Figure 1. A signed random heap with 10 generators. The roof elements are shaded. 

The interpretation of this construction should be clear. Dropping a + (resp. — ) piece in the k th column 
corresponds to adding the generator g^ (resp. g^ ) to the random walk. The fact that gt and gt+i don't 
commute is reflected in the fact that pieces have sticky corners. In section 2, we describe a bijective 
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correspondence between signed heaps with m generators and reduced words in LJ- rn so that a random 
walk on hJ-'m corresponds to the growth procedure of the random heap detailed above. Our bijection differs 
from the one presented in U, which is for the locally free group without periodic boundary conditions. This 
formalism allows us to use the terms random heap process and random walk on hT m interchangeably. 

This geometrical picture of the random walk on hJ- m leads one to consider the roofT(W n ) of the random 
heap, which is the collection of removable pieces. A piece is said to be removable if there are no pieces lying 
above it in its own column or the two adjacent columns. We shall say that a piece x blocks piece y if x must 
be annihilated before y can become part of the roof. The roof is an important feature of the random heap 
because its size controls the expected rate of growth of the heap. Indeed, if the roof has cardinality k, then 
the next piece to fall will annihilate a piece in the heap with probability 5—, and will increase the size of 
the heap with probability 2 ™~ , The long term rate of growth of the heap is thus controlled by the long 
term expected size of the roof. To make this notion precise, we define the drift of the random heap process 
to be £ = linin^oo -^K(#W n ) where 4t=W n denotes the number of pieces in the heap. It follows easily from 
the above observation that £ = 1 — linin^oo EffiZOj^li ; a result given in 0J. 

We prove the following: 

Theorem 1.1. For m > 4, lirn„^oc, ^E(#T(W„)) < 0.32893. 

This result is surprising, since there are heuristic reasons described in 0] for suspecting that the limit should 
converge to | as m — > 00. In particular, one may consider a random walk on the semi-group which has 
generators go, gx,..., g m -i satisfying the relations above, but does not include the inverses of the gt's. This 
random walk W n has an interpretation as an unsigned random heap process. Specifically, we drop pieces 
uniformly over m columns as before, but now we forbid pieces from annihilating. Restricting attention to 
the roof of the unsigned heap, it is easy to see that this process is an irreducible Markov chain and we have 
the following: 

Proposition 1.2. (Proved in ^) If W n is the unsigned random heap process with m generators, then 
lim^oo ±#T(W n ) = § a.s. 

The proof given in 0] may be simplified. If one defines X n> i to be the number of roof elements in the i th 
column after the n th particle is dropped, this process is a Markov chain on {0, 1} with transition probabilities 
Pi : o = 1 — = — and po,i = 1 — Po.o = The stationary distribution for this 2-state chain is (2/3,1/3), 

so lim JWOO E (#T(W n ) \ = lim rwoo Y^iLi ^(X n .i) = y- The claim now follows from the ergodic theorem. 
This elegant proof was discovered by Gabor Pete. From this nice result for unsigned heaps one might suspect 
that the same stationary behavior should be exhibited by the signed process, at least in the limit m — ► 00, 
a conjecture expressed in j^j. 

The paper is organized as follows. First we describe the bijection between elements of LjF m and random 
heaps. Then we prove that for m > 4, 

1 " 

(1) limsup # T (Wk) < 0.32893 a.s. 

n — >oo mn 

h—l 

Finally we shall deduce that if one allows heaps to be infinite, the random heap process has a unique 
stationary distribution v and the finite dimensional distributions of W n converge to the f.d.d.'s of v as 
n — > 00. From this fact, we deduce the same upper bound for the space average, proving the theorem. We 
conclude with some open problems. 

2. The Bijection between Heaps and Words 

We begin by giving a formal definition of a heap, and introducing some notation. 
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Definition 2.1. A heap H is a finite union of pieces (vi,Wi,Oi) with vi G Z + , Wi € Z m and cr^ G {+, — }, 
which satisfy the following conditions: 

1. If (v a ,w a , o~ a ) G H and v a > 1, then there exists (vp,Wp,o~p) £ H such that vp — v a — 1 and 

W/3 G {U>a - l,WQ,W a + 1}- 

2. If (Uq,, w q , a a ), (vp, wp, op) G H satisfy v a — Vp + 1 and w Q = wp, then cr Q = erg. 

3. If (v a , w a , o a ), (vp, wp, up) G H satisfy w a = wp + 1, then w Q ^ wg. 

If if is a heap, the elements of the roof may all be removed to construct a new heap if'. We shall say that 
the roof of H' is the 2 nd level roof of H. Similarly, we may define the 3 rd level roof and so forth. It is easy 
to see that a heap may be determined uniquely by specifying the horizontal positions of the pieces in each 
of its roofs. 

The concept of a roof also makes sense for words representing elements of \jT m . Specifically, if a is an 
instance of gi or g^ 1 in the word W , we say that a is removeable if it commutes with all the letters in W 
occuring to its right, and no other instances of gi or g^ 1 occur to its right. The roof of W consists of the 
collection of removeable letters. Omitting the removeable letters from W to obtain W', we define the roof 
of W to be the 2 nd level roof of W. Higher level roofs are defined analogously. If W = e, the unit in L!F m , 
we say that all the roofs are empty. Observe that two words with the same roofs must represent the same 
element in \jT m , although words with different roofs may also be equal. To obtain a unique representation, 
we introduce the following definition. 

Definition 2.2. A word W representing an element of \tT m is said to be in normal form if the following 
conditions are satisfied. Here a is an instance of gi or g~ l in W , and b is an instance of gj or gj 1 . 

1. If ab = e, then a and b are not members of the same or adjacent roofs. 

2. If a and b are members of the same roof and i < j, then b occurs to the right of a in W . 

3. If a is a member of a higher level roof than b, then b occurs to the right of a in W. 

The fact that every element of L^F m may be represented by a unique normalized word follows directly from 
the commutation relations defining LJ> n . By identifying the respective roofs of a heap and a normalized 
word so that a + (resp. — ) piece in column k corresponds to an instance of the generator g^ (resp. g^ 1 ) one 
obtains the desired bijection, see Figure [5] With this identification, the random growth process of a heap 
described in the introduction corresponds to a random walk on \jT m . 
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Figure 2. The pieces in the second figure arc numbered to indicate their roof membership. The word corresponding to this 

heap is g~ 1 g^ 1 go949ig39295- 

3. Bounding the Time Average 

We now give a proof of inequality l[T)l. Consider a single column of a random heap, and when each piece is 
dropped record a "1" if the column contains an element of the roof and a "0" otherwise. This generates a se- 
quence of 0's and l's, say X\ , X%, A3, . . . , and to prove Q it is enough to show that lim sup,^^ — ^fe=i -^fc — 
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0.32893 a.s. The sequence Xi,X 2 ,X 3 . . . may be expressed more compactly as (0) , (0) 1 , (1) 2 , (0) 2 , . . . 
where denotes the length of the k th sequence of l's and similarly for (0) fc . The random variable (0) is 
distinguished in that its value may be zero. Our approach is to show that 



(2) lim sup — f^ k=1 J^*-i — < 0.32893 a.s. 



ELi (i) fc 

We shall do this by constructing a.s. upper and lower bounds for 



^ oc 1 n 

(3) limsup- (l> fc and hminf - ^ (0) fc , 



n^oo n — ti— >oo n 

k=l fc=l 



respectively. 

An observation fundamental to our proof is that a sequence of 0's may be terminated in two different ways. 

1. A piece may fall in the distinguished column. 

2. All the pieces blocking the highest piece in the distinguished column may be annihilated. 

We shall say that a sequence builds upward if it is terminated via method 1 and otherwise it backtracks. 
Similarly, a sequence of l's may be terminated in two ways. 

1. A piece may fall in a column adjacent to the distinguished column. 

2. All the pieces in the distinguished column blocking the highest piece in an adjacent column may be 
annihilated (or all the pieces in the distinguished column may be annihilated if the adjacent columns 
contain no pieces). 

As before, we say that a sequence of l's builds upward if it is terminated by method 1, and that it backtracks 
otherwise. 

The first task is to construct an a.s. upper bound for limsup,^^ — Efc=i (-Ofcj the following outlines our 
approach. The key idea is based on the observation above. In general, if a sequence of 0's builds upward 
we expect the subsequent sequence of l's to be significantly shorter than if the sequence of 0's backtracks. 
For in the first case the subsequent sequence of l's may be easily terminated via either backtracking or 
building upward. Whereas in the latter case the sequence of l's will likely begin at a time when there is a 
stack of pieces directly below the roof element in the distinguished column so it will be much more difficult 
to terminate the sequence of l's via backtracking. We will describe this phenomenon precisely in Lemma 
13.11 then will show in Lemma 13.31 that sequences usually build upward. Henceforth, we shall say that a 
sequence of l's starts from a short position if it is initiated at a time when there is no piece directly below 
the roof piece in the distinguished column. Otherwise we say that the 1 sequence starts from a long position. 
Also, we shall refer to the distinguished column as column and label the columns to the right 1,2,3, . . . 
modulo m. Unless otherwise stated, we assume without comment that m > 4. 

Lemma 3.1. E((l) fc |(l) fc starts from a short position) = (y2 — l)m and E({l) fe ) < ^ . 



Proof. Consider the following random walk procedure on Z. Start the walk from I > 0, and between 
consecutive steps of the walk wait i.i.d. times tj with distribution P(n > s) = (j^E^-) ■ Also, immediately 
prior to taking each step flip a coin and with probability 2/3 stop the walk. If the walk is not stopped, then 
move 1 unit to the left or right with equal probability. Stop the walk when it reaches the origin, if it has 
not been stopped previously. It is easy to see that this process is equivalent to the process that determines 
the length of (l) fe . The position of the random walk corresponds to the number of consecutive pieces at the 
top of column 0. Setting I — 1 we find: 



(4) 



E((l) fc |(l) fc starts from a short position) = E(r i )E(A r A T), 
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where P(iV > k) = (i) fc and T is the hitting time of for a simple random walk starting from 1. It follows 
from Lemma 3.3 in Chapter 3 of [2] that 

(5) P(T > 2k - 1) = P(T > 2k) 

so, since N and T are independent, 

(6) 



E( W Ar) = i + gO*) (j 



2fe-l 



2// 



= 3(V2-1). 



The above sum is evaluated using the identity YX=o ( kO** = ( X _ 4i)- 1 / 2 which is valid for \t\ < j. Since 
E(n) = f we obtain E((l) fe ) = (%/2 - l)m. From any starting position, E((l) fe ) < E(n)E(N) = f . □ 
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Figure 3. The shaded pieces in each figure are both realizations of the symbol cy^ 1 . 

Next we bound the probability that a sequence of O's backtracks. We say that a sequence of O's starts from 
a short position if it is initiated at a time when there is only one piece blocking the highest piece in column 
zero. Then it is clear that the probability of backtracking is maximized when the series starts from a short 
position. In what follows, we write 

(7) cP , d d , V , % , rf* 3 , n B , V D , °£ , rf^ , rf^ , cr^ , , J 

to denote possible configurations for the pieces blocking the uppermost piece in column zero. In each of 
the above pictures, the lowest piece is assumed to lie in column zero and we agree that exactly those pieces 
which block this piece are shown. The position of a piece indicates its column and which of the other pieces 
it blocks. It does not necessarily indicate its height relative to the others, see Figure [21 Two different 
realizations of a given symbol are said to represent the same configuration of blocking pieces. We shall use 

the notation P( c tP) to denote the probability that a sequence backtracks from the starting position c tP, 
and so forth. To see that this notation makes sense, we appeal to the following lemma. 

Lemma 3.2. The configuration of pieces blocking the highest piece in column zero determines both the 
probability of backtracking and the probability that this piece will be annihilated at any later time. 

Proof. Consider random heaps W 1 and W 2 , and suppose that at time n the pieces blocking the highest piece 
in column zero of W„ and W 2 have the same configuration. Label the highest pieces in column zero x\ and 
X2 respectively. Construct a coupling between W 1 and W 2 as follows. Select a piece uniformly over the 2m 
possible choices and add it to to obtain W^ +1 . Now add a piece in the same column to W 2 , and choose 
its sign according to the following rules. 

1. If the newly added pieces both land on top of roof elements, fix the sign of the second piece so that 
both annihilate or neither annihilate. 

2. Otherwise, choose the sign of the second piece to be the same as the sign of the first piece. 
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Observe that at time n + 1, either both x\ and X2 will have been annihilated or the pieces blocking x\ and 
X2 will again have the same configuration. Iterating this procedure, the proof is complete. □ 

We now estimate P(cP). Observe that by symmetry we have P(cP) = PC-b). 
Lemma 3.3. 0.137457 < P(cP) < 0.14599. 

Proof. Conditioning on the next piece to fall in either column -1,0,1 or 2 we obtain: 

(8) P(rf) = \n V) + \ ■ + \ ■ 1 + ip ( D B) + ip [f) . 

Now expanding the terms on the RHS and using symmetry: 

1 _ / R\ 1 _ / rP\ 1 . 



+ 64" V n 7 ' 32" V n 7 ' 40" 

The lower bound for P(rP) is obtained by setting the unknown P(-) terms on the RHS equal to zero. For 
the upper bound, we use the following lemma. 

Lemma 3.4. For any initial heap, call p the probability that a specific piece in the roof is annihilated at any 
future time. Then i < p < ~. 

Proof. Conditioning on the next piece to fall in any of the three relevant columns, we obtain 

1 1 5 o 

(10) _<,<_ + _,». 

Thus (5p — l)(p — 1) > 0, and we deduce that p < i. □ 
It follows that 

p(^).r(J).r(^),r(/) < i W 

(11) pfl^YrY'tf) < |f(cP). 

Substituting these approximations into ^ yields the upper bound. □ 

The next task is to construct an a.s. upper bound on limsup n ^, 0O ^ J2k=i (l)fe where (0) , (0)i, (1) 2 : ■ ■ • 
are chosen according to the heap process. We construct (1) 1 , (1) 2 , (1) 3 , ... so that (l) k > (l) fc for all k and 
limsupn^oo ^X)fe=i (l)fc nas a sufficiently small a.s. upper bound. The (l) fe 's are constructed as follows. 
Consider the following four random variables: S, S* , S and L. The variable S shall be given the distribution 
of (l) k starting from a short position. We give S* the distribution of S conditioned on the event that the 
first piece to fall in either column -1, or 1 after the 1 sequence is initiated actually falls in either column 
-1 or 1. Also, S has the distribution of S conditioned on the event that the first piece to fall in either 
column -1, or 1 after the 1 sequence is initiated actually falls in column 0. Finally, L has the distribution 
of a (1)^, whose starting position consists of an infinite tower of pieces in column 0. From Lemma 13.11 
E{S) = (V2 - l)m and E(5*) = E{n) = f , moreover E(L) = E(N) = f . Since ±E(S) + fE^*) = E(5), 
we obtain E(S) = (3a/2 - ll/3)m. 

Suppose now that we have an infinite number of i.i.d. copies of S* , S and L. To construct (l) fc , we first build 
a Markov chain on the state space {s* ,§,£}, and then define (l) fc to be a new independent copy of S*, 
S or L if £fc equals s* , s or £ respectively. We shall ignore the condition (1), > (1), for the time being. To 
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construct take £1 = s* with probability 2/3 and £1 = s with probability 1/3. Once has been chosen, 
choose according to the following transition probabilities. Here p = P(cP). 

s* | s* |(1 -p) S * |(1 -p) 

/ / / 

(12) S * - S | S |(1 -p) ~s |(1 -p) 

\ \ \ 

£ £ p I p 

Then ^ has stationary distribution: 

2(1 -p) 1-p 



3-2p ' 3-2p 3-2p 

Since all the copies of S* appearing among the (l) fc 's are i.i.d., and similarly for the copies of S and L, it 
follows from the ergodic theorem that 

1 " - 

(14) - Wfc ^ *>E(S*) + tt s -E(5) + 7T £ E(L). 

fe=i 

Combining equations (|13|l and (|14|) with Lemma T3. 31 we obtain 

1 n ~ 

(15) limsup - Y < 0.41884m a.s. 



n fe= i 



Now, to deduce the same a.s. bound for limsupj^oo — 53fc=i (-0fc> ^ remains only to demonstrate that the 
(l) fc 's can be constructed so that (l) fe > (l) fc . 

Lemma 3.5. We may construct the sequences and (l) fe with the joint densities specified above so that 

(1)^ < (l)fej an d */Cfe = s *j fieri after (l) k is initiated a piece falls in either column 1 or -1 before one falls 
in column 0. 



Proof. The proof is by induction. For the base case, take (1) 1 = To determine £1, consider the first 

piece to fall in either columns -1, or 1 after (1) 1 is initiated. If it falls in either column -1 or 1, then 
£1 = s*, otherwise £1 = s. Now assume that £1, . . . , and (1) 1 , . . . , (1) have been constructed to satisfy 
the stipulations specified above. If = s* we see by the induction hypothesis that (1)„ +1 must start from 

a short position regardless of whether (0) n backtracks. So take (1)„ +1 = (l) n +i m this case. Define as 
before by considering the first piece to fall in either column -1, or 1 after (l) Jl+1 is initiated. 

If £ n = s or £ n = £, then construct a random variable Xn+i independent of (l) l5 . . . , (l) n and £1, ...,£„ so that 
P(Xn+i = 1) = 1 — P and P(xn+i = 0) = p. If H n denotes the configuration of the heap at the time when (0) n 
is initiated, we know from the discussion immediately preceding Lemma l3~3l that P((0) n backtracks \H n ) < P- 
Thus, we may construct Xn+i so that {(0) n backtracks} C {Xn+i = 0}. Now, if Xn+i = 1, take (l) n+1 = 
(l) n+1 . Construct £„ +1 as before by considering the first piece to fall in either column -1, or 1 after (1)„ +1 

is initiated. If £ n = 0, take £n+i = £■ To construct (l) n+1 , start building a random heap process on top 
of an infinite stack of pieces in the column. Drop pieces in the same order as they fell while (l) n+1 was 
being constructed until (l) n , 1 was terminated. At this point, if the column still contains a piece in the 
roof, continue to drop signed pieces uniformly over the m columns, but now independently of the random 
heap process. Take (1)„ +1 to be the length of the sequence of l's that is generated. It is clear that if (l) n+1 
builds upward, it will equal (l) n+1 , and otherwise it will be strictly smaller. □ 

Thus, we have established the following estimate: 
Lemma 3.6. lim sup^^ i Y%=i Wfc ^ 0.41884m a.s. 
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We now calculate an a.s. lower bound for liminfn^oo — X)fc=i The first step is to estimate the condi- 
tional expectation E((0) fc | (0) fc starts from a short position). In what follows, we shall write E( E tP) to denote 
E((0) fe |(0) fe starts from n n n ). 

Lemma 3.7. 0.85453m < E(rf) < 0.86255m. 

Proof. Conditioning on the first piece to fall in either column -1,0,1 or 2 we obtain 

(16) E(cP) = 2 + ±ECtP) + ±E (^) + \E 

and expanding the terms on the RHS again yields: 

<") S^> - ^ - s E (J) - H + (/) + l E K 

Bounding the unknown terms on the RHS by m yields the upper bound. 

For the lower bound, it is useful to first estimate the related quantities P(cP) and E(cP). Here, P(cP) is 
the probability that starting from the indicated configuration, the piece distinguished by the x symbol will 
be annihilated before a piece falls in column 0. Similarly, E(n H ) gives the expected number of steps before 
either the distinguished piece is annihilated or a piece falls in column 0. 

Lemma 3.8. We have the estimates 0.180115 < P(cP) < 0.1806355 and 0.133939 < P ( n n 8 ) < 0.141677. 
For m = 4, P (cf^) = H^) and for m>h we have 0.133939 < P ( C /=P H ) < 0.141677. 

Proof. Conditioning on the next piece to fall in column -1, 0, 1 or 2 we write: 

P(cf) = Jp(W) + i.0 + i.l + ip(^)+ip(jP) 
(18) = Ip(d a ) + i + ip( n B ) 2 + ip(^)p( n B ). 

Expanding P (n ^ by conditioning on the next piece to fall in column 0,1,2 or 3 gives: 

as) 'Kl^'^l+^i'M'ti'f/ 

The probability that the highest piece in n d aD is annihilated at all is p < 1/5, and the probability that it is 
annihilated by the first piece to fall in columns 0, 1, 2, 3 or 4 is 1/10 (actually with 1/10 replaced by 1/8 if 
m = 4, but this extra precision is not important). Thus, we obtain the inequality: 

(20) i.p( D J)<p( nD « D )<lp 

The quantity P^^pP^ satisfies the same inequality. Combining these inequalities with 119|) and 120fl . the 

bounds on P fcF^) follow. Once these inequalities are established, the bounds on P(n a ) follow from 118( 1 . To 

bound P ^ n cP ^ if m > 5, condition on the next piece to fall in column 0, 2, 3 or 4. Then use the analogues 
of l|2U|) to bound the unknown higher order terms. □ 
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The quantity E(n H ) may be readily calculated from P(n H ). Indeed, let T be the number of steps until either 
the distinguished piece is annihilated or a piece falls in column zero. By applying the strong Markov property 
at T, we obtain 



(21) 

Iterating this procedure yields 
(22) E 



E 



=}■ 



l+^d 2 ) E(d a ). 



l + ^d 2 ) + --- + P(d a ) n E(n H ) 



E(n H ) = (l -P(cP)) m 



and letting moowe deduce that 
(23) 



Analogous expressions clearly hold for E Lp j and E ^^pF J . Using (J2HJ) and condit ionings similar to the 
one leading to (|21|l . we are now able to give lower bounds for the unknown terms in (|17|) . 



(24) 
(25) 
(26) 
(27) 
(28) 
(29) 
(30) 



> E 



E (J) - 

E(n n B) > E 

E( n B°) > E 

E( D n nD ) > E 

E(^P) > E 

E (cf fP ) • 



E 



> E 

> E 



> 



1 + P(n H ) + P(d a ) 2 E(n H ) > 0.994106m 
1 + P(n H )l E(n H ) > 0.96737m 



> E( 



(ri^) + P ( j") [l + Hd 3 )] E(n H ) > 



.995377m 



> E \rp") + P E(n H ) > 0.974408m 



(^)e( 



>E( n ^)+P(^) 



E rf 3 + 



5(d") 



> 0.995377m 

> 0.995377m 
> 0.996374m 



The numerical estimates are generated by using the estimates from Lemma |3.8I Substituting these values 
into equation l|17f) now gives the lower bound for E(cP). □ 

It remains to construct an a.s. lower bound for liminf ra ^ 00 ^ Efc=i Wfc- Let (0) be a random variable 
with distribution equal to (0) fc starting from a short position. Since the distribution of (0) fe given Tik is 
stochastically greater than (0), we can construct (0) 1; (0) 2 , (0) 3 , ... to be i.i.d. copies of (0) so that (0) fe < (0) fc 
for all fc. By Lemma |3.7I and the strong law of large numbers, we know that limsup,^,^ — Efe=i (®)k — 
0.85453m a.s. Hence, we arrive at the estimate 

Lemma 3.9. limsup^^^ i Efc=i — 0.85453m a.s. 



It follows that 
(31) 



1 - 

lim sup — } Xk < lim 



sup ■ 



E 



fe=i 



fc=i 



< 



0.41884 

< 0.32893. 



0.41884 + 0.85459 

By symmetry, this result is valid for any column, and equation follows. 
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4. Bounding the Space Average 

We have shown that limsup„_ >00 ^ JX=l # T ( W k) < 0.32893 a.s. Since < #T(W fe ) < m for all fc, it 
follows that 

1 ™ 

(32) limsup ^E(#T(W fc )) < 0.32893. 

n — >oo TTLTL , _.. 

Thus, to prove that limsup^^ ±E(#T(W„)) < 0.32893 it suffices to show that Iim^oo E(#T(W„)) exists. 
We shall prove the much stronger result that if one allows heaps to be infinite, the random heap process has a 
unique stationary distribution v and the finite dimensional distributions of W n converge to the corresponding 
f.d.d.'s of v in total variation distance as n — > oo. 

We begin by introducing some notation. We shall write Tk(W) to denote the configuration of the roofs of W 
above the (fc + l) s * level. Thus, Ti(W) = T(W), T 2 {W) consists of the I s * and 2 nd level roofs, etc. Infinite 
heaps are determined abstractly by specifying the configuration of their roof at each level. The roofs of 
an infinite heap Woo must satisfy the same geometrical relations as those of a finite heap: Tfe(Woo) must 
determine a finite heap satisfying Definition ^, ll for all fc. Finite heaps are a subset of these generalized heaps 
distinguished by the property that only finitely many of their roofs are nonempty. Since the uppermost k 
roofs of a heap are sufficient to determine the uppermost k — 1 roofs after a piece is added, the random heap 
process extends naturally to a Markov operator, V , on the space of generalized heaps. In what follows 
i?^ shall denote the space of possible configurations of the uppermost k roofs, and Ai(Rk) shall denote the 
metric space of probability measures on R k with total variation distance. We work with the sigma field !Fjj 
on Vt H generated by all sets of the form T^ 1 (</>), where <j) € Rk and k is a positive integer. If /x is a measure 
on fi//, it induces a measure fj, o T^ 1 on R^. We claim that V has a unique stationary distribution on the 
measure space {QhtFh), a fact which will follow readily from the following lemma. 

Lemma 4.1. Let W 1 and W 2 be signed heaps starting from random initial configurations Wq and Wq 
with laws and /1q. Then ix„ = V" '/xj and /x^ = V n /J 1 Q, the laws of W„ and W„ respectively, satisfy 
\\Hn° T7 — /x^ o Tr \\ — > as n — > oo uniformly in [Iq and fi^. Here \\ ■ \\ denotes total variation distance. 

Proof. Construct a coupling between W 1 and W 2 as in Lemma l3~2l Now suppose that at some point pieces 
fall successively in columns 0,1,2, ... ,m — 1 and that the pieces falling in column did not annihilate. Once 
this sequence is observed, the roofs of W 1 and W 2 must agree. Call $fc the event that this sequence appears 
k consecutive times. In order for T^(W*) and Ti{W 2 ) to disagree after the event &k+e is observed, at least 
mk of the pieces comprising the $k+£ sequence must be annihilated. By conditioning sequentially on these 
pieces and using Lemma 13.41 this probability is bounded above by (-|) . Now, given e > we may choose 
N,K large enough so that (~) m < ~ and the probability that §k+1 is not observed in the first N steps is 
less than |. Since ||/x* o T^ 1 - [i 2 n o T^ _1 || is bounded by the probability that T e (W^) and T e (W 2 ) differ it 
follows that for all n > N, o Tf 1 - fi 2 n o 1 1 1 < e. □ 

Taking /Xq = /xo and /Xq = \Xj the lemma implies that as n — ► oo, \\fi n o — fi n +j ° T^\\ —> uniformly 
in j. Thus fi n o Tf 1 converges in A4(Re) to a measure V(_, defined by vi{A) — lim^oo /x„ o T^~ 1 (A) for all 
A € Ri, which satisfies vioTgo T^ 1 = for fc < I. It follows from Kolmogorov's extension theorem that 
there is a unique measure v on Qh satisfying v o T^ 1 = for all fc. This measure is the unique stationary 
measure for V. To see that Vv = v, observe that V induces a continuous function V : A4(Rk) — * M(Rk-i). 
Since V{jx n ° T^ 1 ) = /x„+i T^-x an d II Mn ° — ' IIm™ ° ^fc^i — ^fc-ill "~ * as n — > oo, continuity implies 
that V(yk) = Vk-i- Hence, (Vv) ° T^-i = ^fc-i an( i ^ follows that Vv = v. Thus, v is a stationary measure 
and by Lemma l4.1l it is unique. 
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5. Concluding Remarks and Further Problems 

1. We have shown that for a signed random heap, ^ lim^oo W.(#T(W n )) < 0.32893. However, the exact 
value of the limit is unknown for any m > 4. 

2. The random heap process may be generalized as follows. Consider dropping pieces with sticky corners 
uniformly over the m columns, and when a piece lands directly on top of another piece, flip a coin and 
allow the two pieces to annihilate with probability p. When p = 0, this new process is equivalent to 
the unsigned random heap process, and if p = |, we recover the signed random heap process. Taking 
p = 1, we obtain a random walk on the group 

(go, 9i, ■ ■ ■ , 9m-i ■ 9i9j = 9j9i V« ^ j ± 1 mod m, and gf = e Vi) . 

Numerical simulations for small values of m suggest that the expected size of the roof is monotonic 
in p, but proving this remains a challenge. 

Acknowledgements. Many thanks to Yuval Peres for suggesting this problem and for helpful discussions 
and comments. I also thank Bob Hough for helpful discussions, and Gabor Pete for useful comments and 
permission to include his nice proof in the introduction. 
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